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Abstract.  Half-space  problem  of  weak  evaporation  and  condensation  of  a  binary  mixture  of  vapors  is  investigated  on  the 
basis  of  the  BGK-type  Boltzmann  model.  By  a  systematic  asymptotic  analysis,  it  is  shown  that  the  steady  evaporation  and 
condensation  takes  place  only  if  the  parameters  that  characterize  the  state  of  the  condensed  phase  and  that  of  the  vapors  at  a 
far  distance  satisfy  one  relation  for  the  condensation  case  and  three  relations  for  the  evaporation  case.  As  an  application,  the 
resulting  relations  are  used  as  the  boundary  conditions  of  the  Euler  system  in  the  study  of  the  two-surface  problem  of  a  vapor 
mixture  for  small  Knudsen  numbers.  This  system  has  two  branches  of  solutions,  which  causes  two  different  solutions  for  the 
same  physical  situation  in  the  continuum  limit.  This  result  is  discussed  in  connection  with  the  ghost  effect. 


INTRODUCTION 

Half-space  problem  of  evaporation  and  condensation  is  one  of  the  most  fundamental  boundary-value  problems 
in  kinetic  theory  and  has  been  intensively  studied  (see,  for  example,  [1,  2]  and  the  references  therein).  One  of 
the  important  aspects  of  the  problem  is  that  the  information  about  the  conditions  for  the  steady  evaporation  and 
condensation  flows  provides  the  fluid-dynamic  equations  with  the  appropriate  boundary  conditions;  that  is,  it  allows 
us  to  complete  the  macroscopic  description  of  vapor  flows  with  the  phase  change  at  ordinary  pressure. 

According  to  the  studies  so  far,  the  conditions  for  the  steady  flows  depend  qualitatively  on  whether  evaporation  or 
condensation  takes  place  for  a  pure  vapor  and  for  a  mixture  of  a  vapor  and  a  noncondensable  gas.  (There  is  also  a 
qualitative  difference  between  subsonic  and  supersonic  condensations.)  In  the  present  paper,  we  will  first  show  the 
corresponding  results  for  a  binary  mixture  of  vapors  in  the  case  of  weak  evaporation  and  condensation.  By  limiting 
ourselves  to  weak  evaporation  and  condensation,  we  can  take  an  analytical  approach  [3]  to  obtain  the  conditions 
explicitly  for  both  evaporation  and  condensation  cases  and  to  clarify  a  qualitative  difference  between  them.  The 
solution  discussed  here  is  closely  related  to  the  suction  boundary  layer  [4]  in  the  conventional  gasdynamics. 

Next,  we  will  study  the  two-surface  problem  of  a  binary  mixture  of  vapors  for  small  Knudsen  numbers  by  the 
Euler  set  of  equations  with  the  conditions  above  being  the  boundary  conditions.  We  show  that  the  Euler  system  has 
two  branches  of  solutions  in  general,  so  that  two  different  solutions  may  occur  for  the  same  physical  situation  in  the 
continuum  limit.  This  result  is  discussed  in  connection  with  the  ghost  effect  [  1 ,  5]  at  the  end  of  the  paper. 


PROBLEM  AND  FORMULATION 

We  consider  a  binary  mixture  of  vapors,  say  vapor  A  and  B,  in  a  half-space  in  contact  with  their  plane  condensed  phase. 
The  condensed  phase  is  kept  at  a  uniform  temperature  Tw,  and  its  surface  (or  the  interface)  is  located  at  X\  =  0,  where  X, 
is  the  rectangular  coordinate  system.  The  mixture  of  vapors  occupies  the  region  X\  >  0  and  is  in  the  equilibrium  state 
characterized  by  the  pressure  temperature  7k.,  concentration  (molecular  number  fraction)  XA  of  species  A,  and 
flow  velocity  ( w..  0. 0)  at  a  far  distance  from  the  interface.  The  Mach  number  of  the  flow  at  a  far  distance,  which 
is  defined  by  =  Iv’ooKffcTko/wioo)-1/2,  is  supposed  to  be  small.  Here  k  is  the  Boltzmann  constant,  ma  (a  =  A,B) 
is  the  mass  of  a  molecule  of  species  a,  and  ma 0  =  mAXf  +  mHXf,  with  X^  =  1  —  XA.  We  will  investigate  the  steady 
behavior  of  the  vapors  on  the  basis  of  kinetic  theory,  mainly  aiming  at  deriving  the  relation  among  the  parameters 
that  characterize  the  state  of  the  condensed  phase  and  that  of  the  vapors  at  a  far  distance  in  order  that  the  steady 
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evaporation  or  condensation  takes  place.  We  assume  that  (i)  the  behavior  of  the  vapors  is  described  by  the  Garzo- 
Santos-Brey  equation  [6]  (GSB  model)  and  that  (ii)  the  velocity  distribution  of  the  molecules  of  species  a  (a  =  A,  B) 
outgoing  from  the  interface  is  the  half-Maxwellian  at  rest  with  pressure  p":  and  temperature  Tw.  Here  p“  is  the  partial 
pressure  of  species  a  in  the  mixture  saturated  in  contact  with  the  condensed  phase  at  temperature  Tw.  Hereinafter,  we 
call  p“  the  saturation  pressure  of  species  a  and  pw  =  pAv  +  pB  the  total  saturation  pressure,  respectively.  In  the  sequel, 
the  Greek  letters  a  and  /3  are  symbolically  used  to  represent  the  species  of  vapors,  i.e.,  ct.fi  =A,B. 

Let  us  denote  by  c*£  (or  c*£,)  the  molecular  velocity  and  by  (re^c*  3^2)/“  the  velocity  distribution  function  of 
molecules  of  species  a,  where  nfv  =  pA,/kTw  and  c*  =  (2kTw/mA)1/2.  Let  xi  =  X\/{s/n/2)£  be  the  dimensionless 
space  coordinate,  where  £  is  the  mean-free  path  of  a  molecule  in  vapor  A  that  is  in  the  equilibrium  state  at  rest  with 
pressure  p^,  and  temperature  Tw.  Then  the  present  problem  is  written  as  a  half-space  boundary-value  problem  for  fa : 


/ 
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Pwin/m01)  3/2exp(-waC2)  (Ci  >  0),  atxi=0, 
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where  ma  =  ma  /mt,  =  Pw/Pw’  P 2  =  P°°Xff  /  pA„  t„  =  Too/Tw,  and  v«,  =  v^o/c*.  Here  is  a  given  function  of 
£  with  dimensionless  macroscopic  quantities  and  rnH  being  parameters  (see  [6]).  K^a  are  those  constants  related  to 
collision  frequencies  that  are  symmetric  in  a  and  /3  and  are  normalized  as  KAA  =  1 .  In  the  present  paper,  we  denote 
by  nAna,  TwTa,  and  (c*v“,0,0)  the  molecular  number  density,  temperature,  and  flow  velocity  of  species  a,  and  by 
nAn,  TWT ,  and  (c*v,  0. 0)  their  counterpart  of  the  mixture,  whose  definitions  are  omitted  here  (see  [6]).  The  notations 
pApa  for  the  pressure  of  species  a,  pAp  for  the  pressure  of  the  mixture,  Xa  for  the  concentration  of  species  cc,  and 
(c*vi>,0,0)  for  the  average  velocity  of  the  mixture  based  on  the  particle  flux  are  also  used.  They  are  related  to  the 
number  densities,  temperatures,  and  flow  velocities  as  pa  =  hafa,  p  =  fit ,  Xa  =  na /n,  and  w  =  XAvA  +  XBvB . 


ASYMPTOTIC  ANALYSIS  —  SLOWLY  VARYING  SOLUTION 


We  carry  out  an  asymptotic  analysis  of  the  boundary-value  problem  (1)  —  (2b)  in  the  situation  where  the  evaporation  or 
condensation  takes  place  only  weakly.  Following  the  convention,  we  denote  by  e  the  small  parameter  of  the  problem: 
e  =  |voo|.  Since  there  is  no  geometrical  characteristic  length,  the  solution  apparently  has  the  length  scale  of  variation 
of  the  order  of  the  mean-free  path  £  [the  length  scale  of  0(1)  in  the  dimensionless  coordinate  x\  |.  Nevertheless,  we 
consider  a  solution  having  another  length  scale  of  0(£/e)  motivated  by  the  suction  boundary  layer  in  the  classical 
gasdynamics  (see  [3,  4]  and  Chaps.  4  and  5  in  [1]).  We  call  such  a  solution  a  slowly  varying  solution,  following  [1], 

Putting  aside  the  boundary  condition,  we  first  consider  the  solution  of  Eq.  (1)  such  that  d fa /dx\  =  0(/“e)  and 
derive  fluid-dynamic-type  equations  that  describe  its  behavior.  To  do  this,  we  rescale  Eq.  (1)  by  introducing  a  new 
space  coordinate  y  =  xi  £  and  seek  the  solution  in  a  power  series  of  e,  i.e.,  fa  =  /“^  -|-  f".t'  H - .  Correspondingly, 

we  expand  the  macroscopic  quantities  as  h  =  h({)]  +  /q^e  -I - ,  where  h  =  n,nA,nB,v,vA,vB,  etc.  Here  we  assume 

%  =  v“Q)  =  0  because  of  the  physical  situation  under  consideration.  By  substituting  these  expansions  into  the  rescaled 
Eq.  (1)  and  by  equating  the  terms  of  the  same  order  of  e,  we  can  obtain  the  coefficient  functions  (n  =  0.  1,2,...). 
There  are  two  things  to  remark  at  this  stage.  The  first  is  that  is  a  local  equilibrium  distribution  at  rest: 


/(“)  =”“o)(7r^(o)/,»“)  3/2exp(-,»%2/7(o)). 


(3) 


Thus  the  temperature  is  common  to  species  at  the  lowest  order.  The  second  is  that  the  conservation  laws  must  hold  at 
each  order  of  e,  from  which  the  fluid-dynamic-type  equations  are  obtained.  Here  we  show  the  equations  up  to  G(e2). 
The  discussions  in  the  next  section  will  be  done  on  the  basis  of  these  equations. 
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dy  2  w  w  dy  ’ 

2  mBKBA'  P(0)=n(0)T(0),  P(i)  =  «(1)f(0)  +  «(0)7(1) ,  (6) 


and  A  =  A(0)  +  A(j)£  +  •  •  •  and  Ka  =  K +  X“j£  +  •  •  • .  Equations  (4a)  and  (4b)  form  a  closed  set  of  equations  for  Pfoj, 
X^,  7(o),  and  vv(i)-  Equations  (5a)  -  (5c)  form  a  closed  set  of  linear  equations  for  p(\),  X^,  Tjj),  and  m>(2). 

Next  we  consider  the  boundary  condition  at  the  interface  for  Eqs.  (4a)  -  (5c).  At  the  leading  order,  as  is  obvious 
from  Eq.  (3),  can  satisfy  the  kinetic  boundary  condition  (2a)  by  setting  =  1,  pB(y  =  pB,  and  7(0)  =  E  Hence 
the  following  condition  is  obtained  as  the  boundary  condition  for  Eqs.  (4a)  and  (4b): 

P(0)=Pw=  1+Pw>  x^)=X^  =  l/(l+pB),  f(0)  =  l,  at  y  =  0.  (7) 

At  the  higher  orders,  however,  the  slowly  varying  solution  is  seen  not  to  match  Eq.  (2a)  because  of  the  flow  of 
0(e).  Thus  we  introduce  a  correction  to  the  solution  in  the  layer  with  the  thickness  of  a  few  mean-free  paths 
adjacent  to  the  interface.  Such  a  correction,  which  we  denote  by  /“,  has  the  original  length  scale  of  variation,  i.e., 
d /“ /dx\  =  Off").  We  seek  the  correction  fg  in  a  power  series  of  e:  =  f^£  +  f£ 2£2  +  '  •  Substitution  of 

the  expansion  /“  =  /".  +  If",  fx ile  4 —  ■  into  Eqs.  (1)  -  (2b)  gives  a  half-space  boundary-value  problem  of  the 
linearized  GSB  equation  for  fjj!  (n  =  1,2,...);  and  the  analysis  of  the  problem  for  finally  leads  to  the  following 
relation  [7]: 
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where  mw  =  XA  +  mBXB  and  XB  =  1  —  XA.  Here  ay,  aB,  ax,  aA,  aA,  a\,  cy,  cj.  and  cx  are  constants  depending  on  XA 
as  well  as  mB /mA,  KBB ,  and  KAB .  Equations  (8)  -  (10)  are  the  boundary  conditions  at  the  interface  for  Eqs.  (5a)  -  (5c). 


CONDENSATION  AND  EVAPORATION  SOLUTIONS 


We  start  with  observing  the  properties  of  the  solution  of  Eqs.  (4a)  and  (4b)  with  (7)  and  the  condition  at  infinity: 

p  — >  p„  =  pA  +/3b,  XA^XA  =  pA/po o,  f  — >  too,  vv  — +  i£,  as  v  — *■  °°.  (11) 

Here,  in  the  last  condition,  the  plus  sign  is  for  the  evaporation  case  and  the  minus  sign  is  for  the  condensation  case. 

First  it  is  readily  seen  from  the  first  equation  of  (4a)  that  p(0)  is  a  constant.  Thus  from  the  first  conditions  in  Eqs.  (7) 
and  (11)  the  difference  of  /A,  from  pw  is  of  Ole)'. 


Poo-pw  =  0(e).  (12) 

Second,  from  the  second  equation  of  (4a),  i  )  is  also  a  constant.  Thus  the  integration  of  Eq.  (4b)  gives 
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where  A,w  is  defined  by  the  first  equation  of  (6)  with  X“  replaced  by  X" .  Equation  (13)  shows  that  T(0)  and  XA^  are 
monotonic.  Since  Dab ,  A(o),  and  T0j  are  positive,  it  also  shows  that  when  (  >  0,  i.e.,  the  case  of  evaporation, 

Xf())  and  7(0)  are  constants.  Otherwise  they  diverge  at  infinity.  Thus  from  the  second  and  third  conditions  in  Eqs.  (7) 
and  (11)  the  difference  of  X'f  from  XA  and  that  of  t„  from  1  are  also  of  O(e): 

7,-l  =  <9(e),  X*-X*  =  0(e),  for n(0)w(1)  >  0.  (14) 

On  the  contrary,  when  «(o)W(i)  <  0,  i.e.,  the  case  of  condensation,  such  restrictions  do  not  occur.  In  fact,  the  (implicit) 
expressions  of  7)0j  and  Xfy)  are  obtained  for  any  values  of  T,„  and  Xf  by  using  the  monotonicity  (see  Sec.  5.7  in  [1]). 

To  summarize,  there  are  three  restrictions  for  the  evaporation,  while  there  is  only  one  for  the  condensation.  These 
give  those  relations  among  the  parameters  at  the  interface  and  at  a  far  distance  that  must  hold  in  order  that  the  steady 
evaporation  or  condensation  takes  place.  The  relations  are  derived  below  by  the  use  of  Eqs.  (5a)  -  (5c)  with  (8)  -  (1 1). 
The  case  of  condensation  In  this  case,  there  is  only  one  restriction,  which  is  on  the  pressure.  To  see  it  more  closely, 
we  consider  the  equation  and  boundary  condition  for  [the  first  equation  of  (5a)  and  (8)].  Obviously  /i,  i  is  a 
constant,  and  thus  from  the  conditions  (8)  and  (11)  we  have  the  relation 
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In  the  meantime,  dT^/dy  and  dXf^/dv  aty  =  0  can  be  expressed  in  terms  of  the  parameters  at  infinity  if  one  integrates 
Eq.  (4b)  once  and  takes  into  account  the  condition  at  infinity.  Substitution  of  the  resulting  into  Eq.  (15)  gives 
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Equation  (16)  is  the  condition  for  the  steady  condensation. 

The  case  of  evaporation  In  this  case,  there  are  three  restrictions:  the  restrictions  on  the  pressure,  concentration,  and 
temperature  [Eqs.  (12)  and  (14)].  The  relation  for  the  pressure  is  readily  obtained  from  Eq.  (15)  as 


(. poo-pw)e  1 


-avpwW(X) 
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because  X^  and  7(0)  are  constants.  The  relations  for  the  concentration  and  temperature  are  to  be  obtained  from 
Eqs.  (5b)  and  (5c)  with  constant  XA^  and  7j0)  and  conditions  (9)  —  (11).  The  discussion  on  Eqs.  (5b)  and  (5c)  parallel 
to  the  first  half  of  the  present  section  shows  that  and  Xf  ,  are  constants.  Thus,  from  Eqs.  (9)  -  (1 1),  the  following 
relations  for  the  concentration  and  temperature  are  obtained: 


(X^-X^1  =-X^-«v)w(1)  |y=o,  (f.-lJe-^-cv^d)^. 

Equations  (17)  and  (18)  can  be  transformed  into  the  following  relations: 

Poo  Voo  X^  ,  ,  A  .  Voo  Too  Voo 

—  =  1  -av  — ,  — r  —  1  (a\--civ)  — ,  —  =  1— cy  — . 

Pw  -^w  C*  J-w  C* 


(18) 


(19) 


These  are  the  conditions  for  the  steady  evaporation. 


APPLICATION  TO  THE  TWO-SURFACE  PROBLEM 

Consider  a  mixture  of  vapors  A  and  B  between  two  parallel  plane  condensed  phases.  One  of  them  is  kept  at  temperature 
7}  while  the  other  is  kept  at  temperature  7//.  The  interface  of  the  former  is  located  at  Xi  =0  and  that  of  the  latter  is 
at  Xi  =  D.  The  saturation  pressure  of  species  a  at  temperature  7/  (or  7//)  is  pf  (or  p:f) ),  and  the  corresponding  total 
saturation  pressure  is  denoted  by  pj  =  pf  +  pf  (or  pn  =  pfj  +  pf;)-  We  investigate  the  steady  behavior  of  the  vapors 


506 


for  small  Kn.  Here  Kn  is  the  Knudsen  number  defined  by  Kn  =  i/D  with  t  being  the  mean-free-path  of  a  molecule  in 
a  pure  vapor  of  species  A  that  is  in  the  equilibrium  state  at  rest  with  pressure  pf  and  temperature  2}. 

According  to  [5],  when  | pii/pi  —  1|  Kn,  the  overall  behavior  of  the  vapors  is  described  by  the  Euler  set  of 
equations  and  its  boundary  condition  is  determined  by  the  analysis  of  the  half-space  problem.  We  consider  this  situation 
and  further  assume  that  the  difference  of  the  total  saturation  pressures  is  still  small,  i.e.,  Kn  <C  | pn/ pi  —  1|  <C  1.  Then 
the  evaporation  and  condensation  are  so  weak  that  the  boundary  condition  is  given  by  Eqs.  (16)  and  (19).  Here  the 
quantities  with  subscript  °°  are  regarded  as  the  values  that  the  solution  of  the  Euler  set  takes  at  the  interfaces.  Since  the 
solution  of  one-dimensional  Euler  set  is  uniform,  each  of  them  takes  a  common  value  at  both  interfaces.  The  subscript 
w  is  to  be  properly  set  as  I  or  II  depending  on  the  situation  (i.e.,  at  which  interface  the  evaporation  takes  place).  Suppose 
that  the  evaporation  takes  place  at  the  interface  I.  Then  the  condensation  takes  place  at  the  interface  II  (Vo«  >  0).  By 
eliminating  /),„ ,  7 U,  and  XA  from  Eqs.  (16)  and  (19),  one  obtains  the  following  relation: 


V  c{  /  pn 


=  1—1  -a 


ji 


+  (^ 
^Dab 


Tj_ 

Xu  LV  '  c{  )  Tu 
11  (1  —  mB /mA)  a 


■* v 


mu/mA 


XA 

xA 


HGI) 


1/2  1 


(20) 


Here  we  added  sub-  and  superscripts  I  and  II  to  the  quantities  to  indicate  whether  they  are  evaluated  at  the  interface 
I  or  II.  The  corresponding  relation  in  the  case  of  condensation  at  the  interface  I  ( vr„  <  0)  is  obtained  by  the  exchange 
between  sub-  and  superscripts  I  and  II.  Because  of  the  difference  between  the  two  cases,  there  are  two  branches  of 
solutions  for  the  Euler  system  in  general.  We  will  come  back  this  issue  later.  Once  Voo | / is  obtained  from  Eq.  (20), 
Poo,  Too,  and  XA  are  obtained  from  Eq.  (19).  In  this  way,  the  solution  of  the  Euler  system  is  obtained. 

Example  We  give  an  example  of  the  solution  of  the  Euler  set  (Euler  solution,  hereinafter)  obtained  by  the  process 
mentioned  above  .  We  consider  the  case  pi/pf  =  3,  pu  j pf  —  2.8,  pfj/pf  =  1.8,  and  7}// 7}  =  1.  The  parameters  of  the 
molecular  model  are  set  as  mB /mA  =  0.5  and  KBB  =  KAB  =  1.  The  total  saturation  pressure  is  higher  at  the  interface  I 
(Pi  >  pu),  and  the  evaporation  takes  place  at  the  interface  I.  For  the  present  parameter  setting,  Xf  =  1/3  =  0.3333  and 
Xfj  =  1.8/2. 8  =  0.6429,  and  the  coefficients  ay,  ay,  aj,  crv,  aff,  and  aAI  are  given  as  ay  =  1.7139,  ctf/Xf  =  0.7045, 
Cy  =  0.3451,  Oy  =  1 .9069,  a1-/-  =  —0.6484,  and  a%  =  0.8421.  These  are  obtained  numerically  by  solving  the  half¬ 
space  problems  occurring  in  the  analysis  of  the  problem  for  fff  With  these  data,  one  can  solve  Eq.  (20)  and  obtain 
v/c*  =  0.02012.  From  Eq.  (19)  the  other  quantities  are  obtained  as  p/ pf  =  2.8965,  T /7}  =  0.9931,  and  XA  =  0.3306 
(or  pA / pf  =  0.9575).  Here  we  omit  the  subscript  °°.  Figure  1(a)  shows  a  comparison  with  the  solution  of  the  GSB 
equation  (GSB  solution,  hereinafter).  As  Kn  decreases,  a  uniform  region  appears  in  the  GSB  solution  and  extends 
from  the  evaporation  side  (I)  to  the  condensation  side  (II).  The  Euler  solution  predicts  the  values  that  the  GSB  solution 
takes  in  the  uniform  region  correctly.  For  small  Kn,  the  GSB  solution  changes  steeply  near  the  condensation  interface, 
where  the  deviation  from  the  Euler  solution  remains  large.  This  region  is  the  suction  boundary  layer  in  the  classical 
gasdynamics,  which  corresponds  to  the  slowly  varying  solution  discussed  in  this  paper.  If  the  slowly  varying  solution 
is  added  to  the  Euler  solution,  the  prediction  becomes  quite  well  even  in  the  region.  Figure  1(b)  demonstrates  it. 

The  limit  of  pu  — >  pj  The  fact  that  there  are  two  branches  of  solutions  of  the  Euler  system  leads  to  a  possibility 


FIGURE  1.  Profiles  of  the  partial  pressures  pA  and  pB  in  the  case  of  pi/pf  =  3,  pn/Pi  =  2.8,  p/z/pf  =  1.8,  and  7///T)  =  1 
{wB /mA  =  0.5  and  KBB  =  KAB  =  1).  The  dashed  line  and  dot-dashed  line  indicate  the  GSB  solution  for  small  Kn.  The  solid  line 
indicates  the  Euler  solution  in  (a)  and  that  with  a  correction  of  the  slowly  varying  solution  up  to  0(e)  in  (b). 
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of  two  distinct  solutions  for  the  same  physical  situation.  We  will  discuss  it  here.  Let  us  consider  the  behavior  of  the 
mixture  in  the  limit  pu  — >  pj.  Because  of  the  scaling  Kn  <C  | pu/pi  —  1|,  the  continuum  limit  Kn  =  0+  is  considered. 
Note  that  evaporation  is  expected  to  take  place  at  the  interface  I  if  p /  >  pu  and  at  the  interface  II  if  pi  <  pu.  Thus,  the 
solution  for  pu  =  pi~0  is  to  be  obtained  from  Eq.  (20)  while  that  for  pu  =  pi  +  0  is  to  be  obtained  from  Eq.  (20)  with 
the  exchange  between  sub-  and  superscripts  I  and  II.  This  difference  gives  two  distinct  solutions  for  the  same  physical 
situation  when  pn  =  Pi-  Figure  2  shows  an  example.  The  parameters  of  the  molecular  model  are  the  same  as  before. 

The  present  analysis  reveals  that  there  are  two  distinct  solutions 
in  the  continuum  limit  when  pj  =  pu.  However,  there  is  another 
class  of  solutions  for  the  same  physical  situation.  The  Euler  sys¬ 
tem,  on  which  the  above  discussion  is  based,  was  derived  under 
the  assumption  that  Kn  -C  \pn/pi  —  1|.  If  one  considers  another 
scaling  Kn  =  0(\pn/pi  —  1|)  and  takes  the  continuum  limit,  the 
behavior  of  the  vapors  is  described  not  by  the  Euler  system  but  by 
the  convection-diffusion-type  system.  The  details  are  discussed  in 
[5].  According  to  the  reference,  there  are  innumerable  solutions 
of  the  convection-diffusion-type  system  in  the  continuum  limit. 

The  family  of  the  solutions  gives  a  smooth  transition  from  one  to 
the  other  branch  of  solutions  of  the  Euler  system  by  changing  a 
“hidden”  parameter  in  the  convection-diffusion-type  system.  This 
parameter  is  related  to  the  infinitesimal  difference  of  the  total  sat¬ 
uration  pressure  between  the  interfaces.  This  subject  forms  a  new 
category  of  the  ghost  effect. 


0 


0.5 


Xr/D 


FIGURE2.  Partial  pressures  //'  and  //’’  of  the  Eu¬ 
ler  solution  in  the  continuum  limit  in  the  case  of 
Pi/ Pi  =  Pn/pf  =  3,  pfi/pf  =  1-8,  and  Tu/ 7}  =  1. 
The  solid  line  indicates  the  solution  for  p//  =  pi  —  0 
and  the  dashed  line  that  for  p//  =  p/  +  0. 


CONCLUDING  REMARKS 

We  investigated  a  half-space  problem  of  a  binary  mixture  of  vapors  bounded  by  their  plane  condensed  phase  on  the 
basis  of  the  GSB  equation.  Under  the  assumption  of  weak  evaporation  and  condensation,  we  considered  the  slowly 
varying  solution  and  derived  a  fluid-dynamic-type  system,  from  the  latter  of  which  the  relations  among  the  parameters 
that  allow  the  steady  evaporation  and  condensation  were  deduced.  The  relations  were  used  as  the  boundary  conditions 
for  the  Euler  system  in  the  study  of  the  two-surface  problem  of  a  binary  mixture  of  vapors  for  small  Knudsen  numbers. 

In  the  present  work,  the  GSB  model  is  adopted.  Exactly  the  same  result  is  obtained  for  other  BGK-type  models  [8] 
with  a  proper  choice  of  parameters,  because  these  models  give  the  same  transport  coefficients  and  the  same  linearized 
equation  as  the  GSB  model.  The  extension  to  the  Boltzmann  equation  for  Maxwell’s  molecules  is  straightforward. 
However,  it  is  not  so  for  the  other  molecular  models.  The  difference  comes  from  the  fact  that  the  thermal  diffusion 
effect  vanishes  in  the  fluid-dynamic-type  system  exceptionally  for  Maxwell’s  molecules  (and  the  BGK-type  models). 
When  the  effect  is  present,  the  temperature  and  concentration  profiles  are  not  always  monotonic,  so  that  the  discussion 
in  the  fourth  section  does  not  work.  The  extension  to  such  a  case  will  be  treated  in  a  forthcoming  paper. 
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